On the detection of Lorentzian profiles in a power spectrum: A Bayesian approach using ignorance priors

Aims. Deriving accurate frequencies, amplitudes, and mode lifetimes from stochastically driven pulsation is challenging, more so, if one demands that realistic error estimates be given for all model fitting parameters. As has been shown by other authors, the traditional method of fitting Lorentzian profiles to the power spectrum of time-resolved photometric or spectroscopic data via the Maximum Likelihood Estimation (MLE) procedure delivers good approximations for these quantities. We, however, show that a conservative Bayesian approach allows one to treat the detection of modes with minimal assumptions (i.e., about the existence and identity of the modes). Methods. We derive a conservative Bayesian treatment for the probability of Lorentzian profiles being present in a power spectrum and describe an efficient implementation that evaluates the probability density distribution of parameters by using a Markov-Chain Monte Carlo (MCMC) technique. Results. Potentially superior to "best-fit" procedure like MLE, which only provides formal uncertainties, our method samples and approximates the actual probability distributions for all parameters involved. Moreover, it avoids shortcomings that make the MLE treatment susceptible to the built-in assumptions of a model that is fitted to the data. This is especially relevant when analyzing solar-type pulsation in stars other than the Sun where the observations are of lower quality and can be over-interpreted. As an example, we apply our technique to CoRoT observations of the solar-type pulsator HD 49933.Comment: 12 pages, 11 figures, accepted for publication in Astronomy and Astrophysic

The connection between stellar granulation and oscillation as seen by the Kepler mission

The long and almost continuous observations by Kepler show clear evidence of a granulation background signal in a large sample of stars, which is interpreted as the surface manifestation of convection. It has been shown that its characteristic timescale and rms intensity fluctuation scale with the peak frequency (\nu_{max}) of the solar-like oscillations. Various attempts have been made to quantify the observed signal, to determine scaling relations, and to compare them to theoretical predictions. We use a probabilistic method to compare different approaches to extracting the granulation signal. We fit the power density spectra of a large set of Kepler targets, determine the granulation and global oscillation parameter, and quantify scaling relations between them. We establish that a depression in power at about \nu_{max}/2, known from the Sun and a few other main-sequence stars, is also statistically significant in red giants and that a super-Lorentzian function with two components is best suited to reproducing the granulation signal in the broader vicinity of the pulsation power excess. We also establish that the specific choice of the background model can affect the determination of \nu_{max}, introducing systematic uncertainties that can significantly exceed the random uncertainties. We find the characteristic background frequency and amplitude to tightly scale with \nu_{max} for a wide variety of stars, and quantify a mass dependency of the latter. To enable comparison with theoretical predictions, we computed effective timescales and intensity fluctuations and found them to approximately scale as \tau_{eff} \propto g^{-0.85}\,T^{-0.4} and A_{gran} \propto (g^2M)^{-1/4}, respectively. Similarly, the bolometric pulsation amplitude scales approximately as A_{puls} \propto (g^2M)^{-1/3}, which implicitly verifies a separate mass and luminosity dependence of A_{puls}.Comment: 18 pages, 12 figures, accepted for A&

Atmospheric parameters and chemical properties of red giants in the CoRoT asteroseismology fields

A precise characterisation of the red giants in the seismology fields of the CoRoT satellite is a prerequisite for further in-depth seismic modelling. High-resolution FEROS and HARPS spectra were obtained as part of the ground-based follow-up campaigns for 19 targets holding great asteroseismic potential. These data are used to accurately estimate their fundamental parameters and the abundances of 16 chemical species in a self-consistent manner. Some powerful probes of mixing are investigated (the Li and CNO abundances, as well as the carbon isotopic ratio in a few cases). The information provided by the spectroscopic and seismic data is combined to provide more accurate physical parameters and abundances. The stars in our sample follow the general abundance trends as a function of the metallicity observed in stars of the Galactic disk. After an allowance is made for the chemical evolution of the interstellar medium, the observational signature of internal mixing phenomena is revealed through the detection at the stellar surface of the products of the CN cycle. A contamination by NeNa-cycled material in the most massive stars is also discussed. With the asteroseismic constraints, these data will pave the way for a detailed theoretical investigation of the physical processes responsible for the transport of chemical elements in evolved, low- and intermediate-mass stars.Comment: Accepted for publication in A&A, 25 pages, 13 colour figures (revised version after language editing

The nature of p-modes and granulation in HD 49933 observed by CoRoT

Context: Recent observations of HD49933 by the space-photometric mission CoRoT provide photometric evidence of solar type oscillations in a star other than our Sun. The first published reduction, analysis, and interpretation of the CoRoT data yielded a spectrum of p-modes with l = 0, 1, and 2. Aims: We present our own analysis of the CoRoT data in an attempt to compare the detected pulsation modes with eigenfrequencies of models that are consistent with the observed luminosity and surface temperature. Methods: We used the Gruberbauer et al. frequency set derived based on a more conservative Bayesian analysis with ignorance priors and fit models from a dense grid of model spectra. We also introduce a Bayesian approach to searching and quantifying the best model fits to the observed oscillation spectra. Results: We identify 26 frequencies as radial and dipolar modes. Our best fitting model has solar composition and coincides within the error box with the spectroscopically determined position of HD49933 in the H-R diagram. We also show that lower-than-solar Z models have a lower probability of matching the observations than the solar metallicity models. To quantify the effect of the deficiencies in modeling the stellar surface layers in our analysis, we compare adiabatic and nonadiabatic model fits and find that the latter reproduces the observed frequencies better.Comment: accepted to be published in A&A, 9 pages, 5 figure

Modelling a high-mass red giant observed by CoRoT

The G6 giant HR\,2582 (HD\,50890) was observed by CoRoT for approximately 55 days. Mode frequencies are extracted from the observed Fourier spectrum of the light curve. Numerical stellar models are then computed to determine the characteristics of the star (mass, age, etc...) from the comparison with observational constraints. We provide evidence for the presence of solar-like oscillations at low frequency, between 10 and 20\,$\mu$Hz, with a regular spacing of $(1.7\pm0.1)\mu$Hz between consecutive radial orders. Only radial modes are clearly visible. From the models compatible with the observational constraints used here, We find that HR\,2582 (HD\,50890) is a massive star with a mass in the range (3--\,5\,$M_{\odot}$), clearly above the red clump. It oscillates with rather low radial order ($n$ = 5\,--\,12) modes. Its evolutionary stage cannot be determined with precision: the star could be on the ascending red giant branch (hydrogen shell burning) with an age of approximately 155 Myr or in a later phase (helium burning). In order to obtain a reasonable helium amount, the metallicity of the star must be quite subsolar. Our best models are obtained with a mixing length significantly smaller than that obtained for the Sun with the same physical description (except overshoot). The amount of core overshoot during the main-sequence phase is found to be mild, of the order of 0.1\,$H_{\rm p}$.Comment: Accepted in A&

Oscillating red giants in the CoRoT exo-field: Asteroseismic mass and radius determination

Context. Observations and analysis of solar-type oscillations in red-giant stars is an emerging aspect of asteroseismic analysis with a number of open questions yet to be explored. Although stochastic oscillations have previously been detected in red giants from both radial velocity and photometric measurements, those data were either too short or had sampling that was not complete enough to perform a detailed data analysis of the variability. The quality and quantity of photometric data as provided by the CoRoT satellite is necessary to provide a breakthrough in observing p-mode oscillations in red giants. We have analyzed continuous photometric time-series of about 11 400 relatively faint stars obtained in the exofield of CoRoT during the first 150 days long-run campaign from May to October 2007. We find several hundred stars showing a clear power excess in a frequency and amplitude range expected for red-giant pulsators. In this paper we present first results on a sub-sample of these stars. Aims. Knowing reliable fundamental parameters like mass and radius is essential for detailed asteroseismic studies of red-giant stars. As the CoRoT exofield targets are relatively faint (11-16 mag) there are no (or only weak) constraints on the star's location in the H-R diagram. We therefore aim to extract information about such fundamental parameters solely from the available time series. Methods. We model the convective background noise and the power excess hump due to pulsation with a global model fit and deduce reliable estimates for the stellar mass and radius from scaling relations for the frequency of maximum oscillation power and the characteristic frequency separation.Comment: 10 pages, 7 figures, accepted for publication in A&

Atmospheric parameters of 82 red giants in the Kepler field

Context: Accurate fundamental parameters of stars are essential for the asteroseismic analysis of data from the NASA Kepler mission. Aims: We aim at determining accurate atmospheric parameters and the abundance pattern for a sample of 82 red giants that are targets for the Kepler mission. Methods: We have used high-resolution, high signal-to-noise spectra from three different spectrographs. We used the iterative spectral synthesis method VWA to derive the fundamental parameters from carefully selected high-quality iron lines. After determination of the fundamental parameters, abundances of 13 elements were measured using equivalent widths of the spectral lines. Results: We identify discrepancies in log g and [Fe/H], compared to the parameters based on photometric indices in the Kepler Input Catalogue (larger than 2.0 dex for log g and [Fe/H] for individual stars). The Teff found from spectroscopy and photometry shows good agreement within the uncertainties. We find good agreement between the spectroscopic log g and the log g derived from asteroseismology. Also, we see indications of a potential metallicity effect on the stellar oscillations. Conclusions: We have determined the fundamental parameters and element abundances of 82 red giants. The large discrepancies between the spectroscopic log g and [Fe/H] and values in the Kepler Input Catalogue emphasize the need for further detailed spectroscopic follow-up of the Kepler targets in order to produce reliable results from the asteroseismic analysis.Comment: 16 Pages, 12 Figures, accepted for publication in A&

On the detection of Lorentzian profiles in a power spectrum: A Bayesian approach using ignorance priors

Aims. Deriving accurate frequencies, amplitudes, and mode lifetimes from stochastically driven pulsation is challenging, more so, if one demands that realistic error estimates be given for all model fitting parameters. As has been shown by other authors, the traditional method of fitting Lorentzian profiles to the power spectrum of time-resolved photometric or spectroscopic data via the Maximum Likelihood Estimation (MLE) procedure delivers good approximations for these quantities. We, however, show that a conservative Bayesian approach allows one to treat the detection of modes with minimal assumptions (i.e., about the existence and identity of the modes). Methods. We derive a conservative Bayesian treatment for the probability of Lorentzian profiles being present in a power spectrum and describe an efficient implementation that evaluates the probability density distribution of parameters by using a Markov-Chain Monte Carlo (MCMC) technique. Results. Potentially superior to "best-fit" procedure like MLE, which only provides formal uncertainties, our method samples and approximates the actual probability distributions for all parameters involved. Moreover, it avoids shortcomings that make the MLE treatment susceptible to the built-in assumptions of a model that is fitted to the data. This is especially relevant when analyzing solar-type pulsation in stars other than the Sun where the observations are of lower quality and can be over-interpreted. As an example, we apply our technique to CoRoT observations of the solar-type pulsator HD 49933.Comment: 12 pages, 11 figures, accepted for publication in Astronomy and Astrophysic

Observations of tides and circularization in red-giant binaries from Kepler photometry

Binary stars are places of complex stellar interactions. While all binaries are in principle converging towards a state of circularization, many eccentric systems are found even in advanced stellar phases. In this work we discuss the sample of binaries with a red-giant component, discovered from observations of the NASA Kepler space mission. We first discuss which effects and features of tidal interactions are detectable in photometry, spectroscopy and the seismic analysis. In a second step, the sample of binary systems observed with Kepler, is compared to the well studied sample of Verbunt & Phinney (1995, hereafter VP95). We find that this study of circularization of systems hosting evolving red-giant stars with deep convective envelopes is also well applicable to the red-giant binaries in the sample of Kepler stars.Comment: Proceedings paper for the J-P Zahn Symposion, Paris, 6 Pages, 2 Figure

Toward a New Kind of Asteroseismic Grid Fitting

Recent developments in instrumentation (e.g., in particular the Kepler and CoRoT satellites) provide a new opportunity to improve the models of stellar pulsations. Surface layers, rotation, and magnetic fields imprint erratic frequency shifts, trends, and other non-random behavior in the frequency spectra. As our observational uncertainties become smaller, these are increasingly important and difficult to deal with using standard fitting techniques. To improve the models, new ways to compare their predictions with observations need to be conceived. In this paper we present a completely probabilistic (Bayesian) approach to asteroseismic model fitting. It allows for varying degrees of prior mode identification, corrections for the discrete nature of the grid, and most importantly implements a treatment of systematic errors, such as the "surface effects." It removes the need to apply semi- empirical corrections to the observations prior to fitting them to the models and results in a consistent set of probabilities with which the model physics can be probed and compared. As an example, we show a detailed asteroseismic analysis of the Sun. We find a most probable solar age, including a 35 +- 5 million year pre-main sequence phase, of 4.591 billion years, and initial element mass fractions of X_0 = 0.72, Y_0 = 0.264, Z_0 = 0.016, consistent with recent asteroseismic and non-asteroseismic studies.Comment: 15 pages, 5 figures, accepted for publication in The Astrophysical Journal; v2 contains minor changes made in the proofs (updated references & corrected typos